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Abstract 

We consider a nanomechanical analogue of a nonlinear interferometer, consisting of two parallel, flexural nanomechanical 
resonators, each with an intrinsic Duffing nonlinearity and with a switchable beamsplitter-like coupling between them. We 
calculate the precision with which the strength of the nonlinearity can be estimated and show that it scales as 1/n^'^, where n 
is the mean phonon number of the initial state. This result holds even in the presence of dissipation, but assumes the ability 
to make measurements of the quadrature components of the nanoresonators. 
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I. INTRODUCTION 

High-precision measurement is an essential component of any advanced technology, be it classical or quantum. In 
the case of the emerging quantum technologies, however, constraints on our ability to make precise measurements 
are imposed by the Heisenberg uncertainty principle. It is often the case that the measurement objective is simply 
to estimate a single parameter of the Hamiltonian of a system [T]. For example, in atomic clocks the objective is 
to estimate the resonant frequency of a given electronic transition f2]. In the case of an optical interferometer, the 
objective is to estimate an optical phase shift produced by some mechanism of interest, which changes the relative 
path length. As the parameter varies, the dynamics of the system changes, and precise determination of the parameter 
requires that the change in the dynamics of the system be resolvable with sufficient sensitivity to the parameter. 

The precision with which the parameter can be determined depends on the initial state of the system, the nature of 
the Hamiltonian describing the system's evolution, and the measurements to be performed on the system. Most work 
on parameter estimations assumes that system quanta are coupled independently to the parameter, meaning that the 
coupling is quadratic in the field variables, which leads to equations of motion that are linear in the field variables. For 
such linear couplings, the optimal precision in a parameter estimate scales as 1/n, where n is the number of system 
quanta used in the measurement, a scaling known as the Heisenberg limit [3]. Achieving the Heisenberg limit with a 
linear coupling requires using an entangled initial state |1] . If one is restricted to using product states, then a linear 
coupling can only achieve a 1/n^'^ scaling, which is called the shot-noise limit or the standard quantum limit. 

It was recently shown that quantum parameter estimation with scaling better than 1/n could be attained by using 
a coupling to the parameter that is nonlinear in field variables [5]; such super-Heisenberg scalings can be obtained 
even with an initial product state [B]. The use of product states, as opposed to entangled states, circumvents the 
difficulty of creating the entangled states and also makes the scheme considerably more robust against the deleterious 
effects of decoherence. In a related development, a measurement of a phase shift of an optical field [16], using an 
adaptive measurement scheme that requires no entanglement, has achieved a Heisenberg-limited sensitivity in terms 
of number of interactions of photons with the phase shifter, rather than just the number of photons. 

Flexural nanomechanical resonators have an intrinsic Dujjing nonlinearity due to extension on bending |7]; tech- 
nology is progressing towards the point where such resonators can be cooled to near their quantum ground state [H]. 
Thus nanoresonators might provide a system in which parameter estimation beyond the 1/n Heisenberg limit could be 
demonstrated, in this case measurement of the nonlinear coefficient for the Duffing nonlinearity. Other candidates for 
the first experimental demonstration of super-Heisenberg scalings in parameter estimation include the measurement 
of phase in a nonlinear optics setting [51 dni [H], the measurement of a magnetic field in atomic magnetometry \n\ . 
and the measurement of atomic scattering strength in coupled Bose-Einstein condensates [131 [TH [TS] . 
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FIG. 1: Quantum circuit representation of nonlinear nanoresonator interferometer. The input nanoresonator modes experience 
a pulsed beamsplitter-like interaction, evolve according to a nonlinear Hamiltonian, and the beamsplitter-like interaction is 
then pulsed on again. We assume that measurements can be made of either the X or the Y quadrature, of one or both output 
modes (denoted "+" and "— "). Though not shown in the circuit, we have also considered the effect of dissipation accompanying 
the nonlinear evolution. 



Quite apart from any fundamental considerations, the ability to make high-precision measurements of the Duffing 
nonlinearity of a nanomechanical resonator might be of considerable practical interest. The Duffing nonlinearity of 
a nanomechanical resonator is an expression of the applied strain |17j . Nonlinear micro-electromechanical systems 
(MEMS) have already been used to make highly sensitive mechanical strain sensors and accelerometers, with appli- 
cations to engineering and biomedical systems |18j . High-precision measurements of the Duffing nonlinearity would 
also have implications for ultra-sensitive nanomechanical mass and force detection [19 . 

We consider two parallel, flexural nanomechanical resonators, each with an intrinsic Duffing nonlinearity and with a 
switchable, electrostatically-actuated beamsplitter-like couphng between them. The measurement proceeds as follows: 
one nanoresonator is excited into a large-amplitude coherent state, the beamsplitter interaction is pulsed on so that 
the coherent-state excitation is split equally between the two resonators, the nanoresonators evolve independently 



under the nonlinearities (and standard linear dissipation), the beamsplitter interaction is pulsed in the same way 
again, and a homodyne measurement of the nanoresonator quadratures is performed. As depicted in Figure [l] this 
scheme effectively realizes a nonlinear interferometer [5D]. We have calculated the precision with which such a scheme 
can estimate the nonlinear coefficient of one nanoresonator and demonstrated that the precision scales as 1/n^/^, 
where n is the mean phonon number of the initial coherent state. 

II. SYSTEM PROPERTIES AND HAMILTONIAN 

Each nanoresonator can be regarded as a thin bar of length / and lateral width a. We are interested in the 
fundamental mode of vibration of each nanoresonator in the lateral direction. Each fundamental mode is described 
by a position- momentum pair, Xi-pi, where i = a,b labels the resonators. With a time-dependent capacitive coupling 
dependent on the displacements from the equilibrium positions, C{xa,Xb), and Duffing nonlinearities characterized 
by coefficients Xii i = ci,b, the system can be described classically by a Hamiltonian 

H,i = -mu^xl + ^ + ^ruLo^xl + |^ + P{t)-C{xa,xi,)V^ + -XamLu^x^ + -Xh^u'^xl , (1) 

where P{t) specifies the coupling voltage pulses. 

In current experiments, the nonlinear coefficient % of a nanoresonator is estimated [71 by measuring the critical 
amplitude Oc at which the forced oscillations become bistable [21 and the quality factor Q of the oscillator and then 
using the relation x = 2\/3/9a^Q. Achievable values of these parameters are Oc — 0.7 nm and Q — 20 000, giving a 
nonlinear coefficient x = 4 x lO^'^ m~^. We use this value of x as a typical value in the following. 

We assume that the nanoresonators are capacitively coupled to nearby bias conducting surfaces in such a way that 
for small displacements, the capacitance can be expanded as 

C(x.,..) = Co(l+^^i±4±'^ + ...). (2) 

Here Cq is the capacitance when the oscillators are at their equilibrium positions. The capacitive coupling must be 
balanced so that there is no net force on the resonators when the coupling is switched on (i.e., no linear terms in 
the expansion). This leaves the quadratic terms as the dominant effect of the coupling. In the quadratic terms, 
d ~ 100 nm is a characteristic lateral separation between the resonators and the other conducting surfaces and / is a 
factor of order unity. Both d and / depend on the specific design of the capacitive coupling. Provided 

CoV^/2mLjd^ = K < cj , (3) 

we can neglect the renormalization of the resonator frequencies during the pulsing of the capacitive coupling, retaining 
only the coupling between the resonators. With these assumptions, the capacitive term in the Hamiltonian (fTj) can 
be replaced by P{t)CQVQ XaXf, / d"^ , which gives rise to the desired beamsplitter coupling. The parameter k, introduced 
in Eq. ([3|, characterizes the strength of the beamsplitter coupling. 
Now we quantize by introducing the operators 



ia = {a + a^)^^, p, = -^(a-a^)^^, (4) 

satisfying the usual commutation relations, and similarly for nanoresonator b. Transforming to an interaction picture 
and using the rotating-wave approximation, we find 

H = hj{a^a)^ + h(3{bhf + hKP{t){o)b + ab^) , (5) 

where 

7 = ^-uxai^xf , (3 = ^^uoxbi^xf , (6) 



with Ax — ^Jhj^/mu) being the half-width of the ground-state wave function. The use of the rotating-wave approxi- 
mation requires that the time scale bt over which the coupling is varied in P(t) be significantly longer than the periods 
of the nanomcchanical resonators. Since 8t ~ k~^, this is the requirement, already introduced in Eq. (l3|, that k <^ w. 
At low temperatures, nanoresonator damping is thought to be mostly due to coupling to a bath of two-level 
systems, but this mechanism is not fully understood. Therefore we do not try to model this mechanism, but rather 



treat dissipation using a quantum optics master equation (with a zero-temperature bath) , with the expectation that 
this generic model of damping provides a reasonable account of the effect of dissipation on parameter estimation. 
Then the evolution of the density matrix describing the state of the two nanoresonators is given by 

p(t) = ~ 1 [iJ, p] + ^ {2apa) - a) ap - pa) a) + ^ (2&p6^ -b^bp- pb"^ b) . (7) 

Experimentally reasonable values for the system properties are I = 2/im, a — 40 nm, m — 10^^^ kg, lo — 2-11 x 
15 MHz = 9.4 X lO'^rad/s, d = 120 nm, Q = 20 000, x = 4 x W^ vcr"^ , Cq = 10 aF, and Vq = 1 V. These correspond 
to Aa; = 240fm, Ti — 4700s~^, 7,/3 = 1.6 x 10~^s^^ (we use lO^'^s^^ as a typical value in the following), and 
K = 3.7 X 10^ s~^. The quantity l/2x(Ax)^ = 5 x 10^^ is roughly the number of phonons required to make the quartic 
nonlinearity as large as the harmonic potential; it corresponds to an oscillation amplitude 1/y^ = 200 nm. 

In Sec. |IV| we analyze estimation of the nonlinear coefficient 7 of oscillator a, assuming that oscillator b has no 
nonlinearity (/3 = 0). Other operating conditions are possible and yield similar results, but the figures in the remainder 
of the paper refer to the (3 = case. We consider a fiducial evolution time t = 10~^ s, so that Tit = 4.7, meaning that 
the effects of dissipation are large, but not overwhelming, and we consider a fiducial initial phonon number n = 10^, 
so that the nonlinear phase shift nji is about 1 rad. We investigate values within about an order of magnitude of 
these fiducial values. Notice that a phonon number n = 10^ corresponds to an oscillation amplitude Axy/2n = Inm. 
This amplitude is close to the value we assumed for Cc, not by accident, but because the two oscillation amplitudes 
quantify, one for free oscillations and one for forced oscillations, the same measure of the relative strengths of the 
nonlinearity and the damping. 

Our assumptions about the switchable beamsplitter coupling require that k <^ us, so that we can make the rotating- 
wave approximation, and that t ^ k~^ = 2.7 /xs, so that we can regard the beamsplitter pulses as essentially instan- 
taneous. Both of these inequahties are well satisfied by the above values. 

III. SYSTEM EVOLUTION IN THE Q REPRESENTATION 

Now suppose that nanoresonator a is excited into a coherent state with amplitude ao, assumed real. The (product) 
state of the two nanoresonators is then given, in the Q representation [22], by Q{a,a*;t = 0) = e^'^^""' /n and 
QiP, (3*;t = 0) = e^l'^l'/TT. If we pulse on the beamsplitter interaction for a time 5t = tt/Ak, the state of the system 
is a product of two equal-amplitude coherent states; in other words, we have a "balanced" beamsplitter. The pulse 
time 6t is sufficiently short that the effect of nonlinearities and dissipation are negligible. We then have 

Q(a,a*;t = (5i)= -e-l"-"°/^l\ g(/3,/3*;i = <5t) = - e-l'^-""/^!'. (8) 

TT vr 

Setting P{t) = in Eq. (Is]) for the time between the beamsplitter pulses, the master equation (|7| can be converted 
into a Fokker-Planck equation for the Q function and solved to give |13] 
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(9) 
(10) 



5a = 2-i{p-q)/Ta. 6b^2p{p-q)/Ti,, (11) 

fa{t) = exp[-rat - 2i-ft{p - q)] , fb (t) = exp[-rbf - 2i(3t{p - q)] . (12) 

The beamsplitter interaction is pulsed on again for time St, giving output quadratures 

a + at±5±6t i{a-a^±bTb^) ,^ „, 

""^^ 71 ' ^ = 71 ■ ^ ^ 



We can calculate the first and second moments of these quadratures using Eqs. ^ and (10 1, 
{X±) = V2 I d^ane{a)Q{a,a*;t)±V2 I d^(3'Jle{l3)Q{(3,(3*;t) , 
(Xl) = -1 + 2 f d^a[JZe{a)fQ{a, a* it) +2 f d^(3[Jle{f3)fQ{f3,f3* it) 
±4 f d^aTZe{a)Q{a,a*;t) f d^(3TZe{(3)Q{(3, (3*;t) . 



(14) 



(15) 



Corresponding moments of the conjugate quadratures Y± are given by the same expressions with the replacement 
TZe{a) -^ Jm{a). 

The evaluation of these moments reduces to the calculation of two integrals, 



V2 fd^aaQ{a,a*;t) = 7^e-(r.*+»c.)/2g.ht+„D./2) _ 
2 I d^a[C{a)]^Q{a,a*-t) = 1 + ^e"^"* ± ^e-'^"*-"^^/2^os(47i + nD4/2) , 
where n = a^, C{a) ~ TZe{a) for the upper sign and C{a) = Tm{a) for the lower sign, and 
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FIG. 2: Magnitude of expectation values of output quadratures X+ and Y+ as functions of the nonlinearity 7 (expressed in 
terms of the nonlinear phase shift njt) and nanoresonator damping r(= Fa = F;,), for the choices n — 10^, (3 — 0, and 
t = 10""^ s. The rapidly oscillating fringes with respect to 7 give rise to the enhanced sensitivity of this parameter estimation 
scheme. Damping leads to a shift in the location of the fringes and to a decay of the quadrature expectations. 



The general results are complicated and thus are not quoted here. In the case of no damping (Ta = = Fh), we 
can show that 



(X±) = y/iie-"''^^-''°'''^-''^/^cos(jt+-sm2jt]±y/7ie- 



is(f3t- 



n(l-cos2/3t)/2cos(/?i+-sin2/3i) , 



{Y±) = V^< 



,-n(l-cos27t)/2 ^ 
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FIG. 3: Uncertainty in the output quadratures X+ and Y+ as a function of the nonhnearity 7 (expressed in terms of the 
nonlinear phase shift wyt) and nanoresonator damping r(= Va ~ Tb), for the choices n = 10^, P ~ Q, and t = 10"'^ s. Squeezing 
and anti-squeezing are observed for low damping, corresponding to shearing of the contours of the Q function in phase space. 
Dissipation suppresses this effect, and coherent-state variances, corresponding to decay to the vacuum, are seen at high damping 
rates. 



and 



{xD 



l+„+" -"(l-cos47t)/2^Qg 



(n 



[A-ft + ^ sin47i) + ^e-"(i-™^4/3t)/2 ^osi^^pt 
^„g-n(2-cos 27t-co. 2/5t)/2 ^^^ (^^^ " /?* + ^ sm 27^ - ^ sin 2/3t) 
±„g-n(2-cos 27t-co. 2/5t)/2 ^^^ |^^^ + /3i + ^ sin 27* + ^ sin 2/3t) , 

^„g-n(2-cos 27t-cos 2/5t)/2 ^^^ |^^^ _ /Ji + !^ sin 27* - " sin 2/3t) 
Tne-"(2-c°>^ 27t-cos 2/3*)/2 ^Qg (^^^ + (it+- sin 27* + - sin 2/3t) . 



■ sin 4/3t 



(22) 



— sm 
2 



in4/3i) 



(23) 



Note that these results could have been calculated directly in the Heisenberg picture with initial coherent states in 
each mode. Making, in addition, the short-time approximation, n{jt)'^ , n(f3t)'^ <^ 1, which still allows nonlinear phase 
shifts 717* Etnd n/3t much larger than unity, we find the same expectation values as for the analogous classical nonlinear 
interferometer: 



{X±) —r ^/n COS njt± ^/n cos n(3t , {Y±) ^ \/n sin njt ± ^/n sin n^t . 

We can define a regime of strong damping by the conditions 
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The conditions on the evolution time are of little consequence because long before they are violated, the oscillators will 
have damped to the ground state. Notice that these conditions allow the case of most interest to us, i.e., Fb = Fa ~ ^n, 



with n'yt ~ Fat allowed to be considerably larger than unity. In the strong damping regime, we have 
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(1±)2, which means that the quadrature uncertainties 
are at the coherent-state level, i.e., AX± = AKj- — 1. 

The dependence of the expectation values of the "+" quadrature amplitudes on the nonlinearity and damping 
rate is shown in Figure |2] similar behavior is displayed by the "— " quadratures. The nonlinearity gives rise to 
rapidly oscillating fringes, and this is the key, as we see in the next section, to the enhanced sensitivity of the 
nonlinear interferometer. Dissipation leads to a reduction in the expectation values and also to a reduction in the 
fringe frequency. Figure [3] shows the uncertainties in X^ and y+, also as functions of the nonlinearity and damping 
rate. In the case of low damping, as the nonlinear phase shift increases, the variances oscillate, corresponding to the 
shearing apart and partial recurrence of the contours of the Q function in phase space. Squeezing and anti-squeezing 
are apparent in the variance oscillations |21], and the squeezing can be quite substantial. Damping suppresses the 
shearing — and, hence the squeezing and anti-squeezing — and suppresses the quantum interference effects that give rise 
to partial revivals. Strong damping leads to a decay of the expectation values and to quadrature variances that take 
on the coherent-state value. At reasonable damping rates, we cannot derive any benefit from the squeezing, though 
neither are we adversely affected by the anti-squeezing of the quadratures. The deleterious effect of damping is almost 
entirely a consequence of reducing the signal carried by the expectation values, not of changing the variances. 



IV. PRECISION OF PARAMETER ESTIMATE 



To analyze the precision of estimating the strength of the Duffing nonlinearity, we specialize in this section to the 
case where oscillator h has no nonlinearity (/3 = 0). We are thus estimating the nonlinear coefficient 7 of oscillator a. 
Other operating conditions are possible, but we focus on this one as a representative possibility in this section. 

We phrase our results in terms of the precision in estimating the related dimensionless parameter 7^, with t regarded 
as fixed. The uncertainty in an estimate of ^t based on multiple measurements of a quantity Z — in our case, Z is one 
of the output quadratures — can be calculated from 



(5(7i) =td-/ = t 



AZ 



AZ 



\d{Z)/dj\ \d{Z)/d{jt)\ 



(28) 



where AZ is the uncertainty in Z. In the case of no damping and again making the short-time approximation, the 
quadrature variances all take on coherent-state values, i.e., AX±,AY± -^ 1. The precision of the estimate of jt thus 
becomes 
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These sensitivities oscillate with the fringes produced by the nonlinear phase shift wyt, but they all have the same 
basic scaling of 1/rt'/'^ with phonon number. This scaling beats the 1/n scaling achievable with a linear Hamiltonian 
and is consistent with the general result fBl for nonlinear Hamiltonians and initial product states. The factor of 
n enhancement compared with the standard quantum limit for linear Hamiltonians is a consequence of the rapidly 
oscillating fringes in the expectation values of the output quadratures. 

From an experimental perspective, the strong damping regime (25 1 is most relevant. In this regime, the quadrature 



variances have coherent-state values, and the derivatives of the expectation values lead to sensitivities 
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The improved l/ii?''^ sensitivity scaling survives in the presence of dissipation, but the absolute sensitivity is degraded, 
and the fringes become more widely separated. For feasible damping rates, the sensitivity is worsened by less than an 
order of magnitude, but if the damping is further increased, the sensitivity diverges, reffecting the absence of signal 
in the quadrature expectation values. 

Figure [4] shows the measurement precision for measurements of the X_|_ and 1+ quadratures as a function of the 
nonlinearity 7 and for two values of the damping rate F = Fq = F;,. Fringe boundaries are located at 777^ = 7?77r/2; 
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FIG. 4: Precision S{'yt) for measurements of the X+ and Y+ quadratures as a function of the nonlinearity 7, expressed as the 
nonhnear phase shift wyt, for the choices n = lO'^, f3 = 0, t = 10"'^ s, and Fa = Fb = F. Zero damping and moderate damping 
cases are shown for each quadrature. For zero damping, fringe boundaries are located at n^t — rmv/2, with the fringes based 
on measurements of conjugate quadratures displaced by 7r/2. Dissipation leads to an overall reduction in sensitivity, and the 
fringes become more widely spaced. 



those based on measurement of the X^ and Y^ quadratures are displaced by 7r/2. As the damping rate increases, 
the overall sensitivity worsens, and the fringes become more widely spaced. These effects can be traced back to the 
reduced-amplitude and reduced-frequency oscillations of the quadrature expectations as a function of the nonlinear 
phase shift. 

The scaling of the measurement precision as a function of n is plotted in Figure |5] Here n is chosen so that njt < 1. 
The precision associated with measurement of the 1+ quadrature is then near its optimal value, away from its first 
fringe boundary at n'yt — 7r/2, whereas the precision associated with measurement of the X_|_ quadrature decreases 
rapidly as it falls from the very poor sensitivity near its central fringe boundary at n"ft — 0. From the log-log plot, we 
can calculate that 6x^ oc vT^I"^ and 8y^ oc n~'^/^, though the extra n~^ in the n~^^^ scaling is due to the sensitivity 
falling from the central fringe boundary, and the true scaling of the optimal sensitivity achievable is n^'^/^. The scaling 
behavior is maintained in the presence of feasible levels of dissipation, although there is a marked deterioration in 
sensitivity. 

In practice, the initial state of the excited nanoresonator would be better described by a displaced thermal state. 
Provided that the thermal width of this state is small compared with the amplitude of the displacement, a condition 
that would be well satisfied by an excitation at the level we are considering, the primary effect of an initial thermal 
distribution would be simply to increase the output quadrature variances, leading to a reduction in sensitivity, but not 
to change in the scaling behavior. In the strong damping regime, contact with a finite temperature bath would result 
in thermal rather than coherent state variances in the output quadratures, again leading to a reduction in sensitivity, 
but not to a change in the scaling with n. Note that the quantum optics master equation for the nonlinearities we 
consider and with a finite temperature bath has been solved analytically using the Q representation [2S] . We conclude 
that the measurement scheme we describe is reasonably robust to increases in temperature. A more severe difficulty 
lies in performing quantum-limited homodyne detection of the output quadratures, though it is conceivable that such 
measurements could be performed using coupled microwave cavities ^26; . 
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FIG. 5: Precision 5["ft) for measurements of the X+ and Y+ quadratures as a function of mean phonon rmmber n in the initial 
coherent state, for the choices 7 — 10~* s~^, /3 = 0, f = 10""^ s, Fa = Fb = F. Plots for three values of the damping constant F 
are shown for each quadrature. These plots correspond to the regime n^t < 1. From the log-log plots, we see that 5x^ oc n~^' ^ 
and 5yj^ oc n~^". The extra n~^ factor for measurement of the X+ quadrature is due to the precision improving as one moves 
away from the very poor sensitivity near the central fringe boundary. 



V. CONCLUSIONS 

We have calculated the precision with which the nonlinearity of a nanomechanical resonator can be estimated, using 
a nanomechanical analogue of a nonlinear interferometer. For an input coherent state, the precision scales as 1/n^", 
a scaling beyond that achievable with a linear coupling even when entangled input states are employed. This scaling 
behavior is maintained in the presence of dissipation, which we modeled using a quantum optics master equation, 
and it is expected that this scheme is reasonably robust to increases in temperature. Quantum-limited homodyne 
detection of the nanoresonator quadratures is, however, a very challenging experimental task. 

An alternative scheme would use a single nonlinear nanomechanical resonator coupled to the field in a supercon- 
ducting microwave cavity, which would act as the second "arm" of an interferometer. The nonlinearity would then 
only be in the mechanical arm of the interferometer. The initial state would be excited by driving the cavity, and 
read-out would be performed by quantum-limited homodyne detection of the cavity output. A beamsplitter-like 
interaction between the cavity and nanoresonator could be realized by driving the cavity on its red sideband [57]. 
This beamsplitter coupling would be continuous, so the analytical results obtained here are not directly applicable. 
Investigating the achievable sensitivity of this cavity-nanoresonator scheme is the subject of ongoing work. 
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